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We present the first lattice QCD determination of the At — >■ A transition form factors that 
govern the rare baryonic decays A;, — >■ A£'^£~ at leading order in heavy-quark effective theory. Our 
calculations are performed with 2+1 flavors of domain- wall fermions, at two lattice spacings and with 
pion masses down to 227 MeV. Three-point functions with a wide range of source-sink separations 
are used to extract the ground-state contributions. The form factors are extrapolated to the physical 
values of the light-quark masses and to the continuum limit. We use our results to calculate the 
differential branching fractions for At, — >■ A£^£~ with £ = e, jj.,T within the standard model. We find 
agreement with a recent CDF measurement of the Ab — >■ A jj.^ fJ,'^ differential branching fraction. 

PACS numbers: 12.38. Gc, 12.39.Fe, 12.39.Hg, 14.20.Mr 



I. INTRODUCTION 



The flavor-changing neutral-current transitions b s ^ and b s can provide strong constraints on models 
of new physics. The effective weak Hamiltonian describing decays of this type has the form ^nSi 
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where the operators that directly give "short-distance" contributions to these decay amplitudes are 
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and the C^ ^ are Wilson coefficients. In the standard model, the Wilson coefficients of the scalar and pseudoscalar 
operators Os,p, as well as those of the opposite-chirality operators O'^, are highly suppressed [S]. Experimental 
measurements of 5 — > s decay rates, angular distributions, and related observables provide constraints on various 
functions of the Wilson coefficients Ci and [9]. In this way, these measurements restrict models of new physics 
and their allowed parameters. Most of the existing studies have focused on mesonic decays such as i3 — ;> fC*7 and 
B — > K^*'>£^£^, for which experiments have reached a high level of precision. To get the most complete set of 
constraints on new physics, it is important to consider many different observables. To this end it is useful to analyze 
also baryonic 6 — >■ s decays such as Af, — A7 and A;, — >■ A£~^£~. The decay Ab — >■ A/i+/j,~ has recently been observed 
using the Tevatron pi) , and is being measured at LHCb. 

One important aspect that distinguishes A^ decays from B meson decays is the spin of the A^ baryon, which in 
principle provides an additional handle on the fundamental interactions. When produced through Z bosons at e'^e" 
colliders, b quarks have a strong longitudinal polarization, and the Af, baryons keep most of that polarization |11H16) . 
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At the Tevatron and the LHC, the Af, baryons produced in proton-(anti)proton coUisions are expected to have some 
degree of transverse polarization |17H19j , which can be measured accurately using the method proposed in Ref . |20j . 
As first mentioned in Ref. [5T] and later studied in detail in Ref. [22], the A(, polarization can be exploited to test 
the "helicity structure" of HcS, that is, to disentangle the contributions from the Wilson coefficients Ci and C'^. In 
practice, this entails measuring an asymmetry in the angular distribution between the At, spin and the momentum of 
a particle in the final state |22j . Even for unpolarized A^ baryons, the spin of the final-state A baryon can also be 
exploited to test the helicity structure of Hoff , as discussed for Af, — A 7 in Refs. p2H?7] and for A^ — A in 
Refs. f25, 26, 28 30 . To this end, an angular analysis needs to be performed for the secondary weak decay A — p7r~. 
Lepton asymmetries for Ab A £^£^ have also been considered 

In order to use these Ab decays to search for new physics, the matrix elements (A7|'Hcff |Ab) or {Ai'^ £~ \7ics\-^b) 
must be determined. For the operators in Eq. ([2]), this then requires the computation of the hadronic matrix elements 
(A| sFfejAb), which are expressed in terms of ten QCD form factors. When using heavy-quark effective theory (HQET) 
for the b quark, the number of independent Af, — > A form factors reduces to 2 JHS'JSS]. Furthermore, in the limit of large 
recoil, soft-coUinear effective theory (SCET) predicts that only one form factor remains ,27, ,36, i37j . The Ab — >■ A form 
factors have been estimated using various models or approximations, including quark models [55tilT| . perturbative 
QCD and sum rules [25l |26l [36l HSI , but have not been determined from QCD previously. In the charm sector, 
some information on the A^, — )■ A form factors is available from the experimental measurement of the semileptonic 
Ac Ae+i^e decay [311l3Sj, and this information has been used to constrain the Ab — > A form factors in Refs. p5ll38j . 
In summary, several estimates of Ab — > A form factors exist in the literature, but a considerable uncertainty remains, 
especially in the low-recoil region where SCET and light-cone sum rules are not applicable. Clearly, first-principles, 
nonperturbative QCD calculations of the form factors are needed, and the method for performing such calculations 
is lattice QCD. 

In this paper, we report on the first lattice QCD calculation of Ab — t- A form factors (we presented preliminary 
results of this work in Ref. [46 ). We use HQET for the b quark, and compute the two form factors that appear. Their 
definitions are given in Sec.|llj Treating the b quarks with HQET on the lattice [47] also leads to several other technical 
simplifications that make the calculation feasible, as will become clear in Sec. [HI A[ For the up, down, and strange 
quarks, we use a domain-wall fermion action [48 50 , which is computationally expensive but provides chiral symmetry 
even at non-zero lattice spacing. Our calculations make use of gauge field ensembles generated by the RBC/UKQCD 
collaborations [51 . These ensembles include 2-1-1 fiavors of dynamical sea quarks, and the lattice parameters used in 
our study are given in Sec. |III B[ We use two different lattice spacings and several different values of the light quark 
masses, which allows us to perform simultaneous extrapolations of the form factors to the continuum limit and to 



the physical values of the quark masses. The data analysis involves several stages, which are explained in Sees. HI C 



|IIID[ |III E[ and [HI F] and estimates of the systematic uncertainties in the form factors are given in Sec. |III G| As a 
first application of our form factor results, in Sec. |IV| we then calculate the differential branching fractions for the 
decays Ab — )• A £~^£~ with £ = e, r in the standard model. The differential branching fraction for Ab — >■ A fX^ 
can be compared to the existing Tevatron data and is of immediate interest for LHCb. Further phenomenological 
applications of the form factors that we have determined are left for future work. 



II. DEFINITION OF FORM FACTORS 

In QCD, using Lorentz symmetry and the discrete C, P, T symmetries, one can show that the matrix elements 
(A(p', s') \ sTb |Ab(p, s)) with F = 7^, 7/^75, q'^a^^, (fcf^^^^ (where q = p — p') are parametrized by ten independent 
form factors (see for example Ref. [5S]). Heavy-quark symmetry, which becomes exact in the limit Wb — >■ 00, reduces 
the number of independent form factors to 2 [55H55] . In the following, when working with HQET, we denote the 
heavy quark by Q. The Ag — ^ A matrix element with an arbitrary Dirac matrix F is then given by |33H35j 

(A(p', s')\ sTQ \Aq{v, 0, s)) = uip\ s') [Fi(p' ■ v) + ^ F^ip' ■ v)]TU{v, s), (3) 

where v is the four- velocity of the Ag, and the two form factors Fi and F2 can be expressed as functions solely of 
p' ■ V, the energy of the A baryon in the Ag rest frame. Here we use the following normalization of states and spinors: 

(A(p,s)|A(p',s')) - 2EA{2nfS,s'6Hp-p'), (4) 
{Aq{v, k, s)\Aq{v, k', s')) = 2i;0(27r)35,,,^3(k - k'), (5) 
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For most of the analysis in this paper, it is convenient to work with the hnear combinations 

F+ = Fi + F2; 
F- = Fi — F2, 



(8) 



instead of Fi and F2. Note that in the hmit mt, 00, five of the ten hehcity-based A;, — ?> A form factors introduced 
in Ref. [36] become equal to F^, while the other five become equal to F„ [see Eq. (2.10) of Ref. [36], where Fi is 
denoted by A, and F2 is denoted by B] . 



III. LATTICE CALCULATION 



A. Two-point and three-point functions 



For the lattice calculation, we work in the Aq rest frame, so that v — (1, 0, 0, 0). The heavy quark, Q, is implemented 
with the Eichten-Hill lattice HQET action [37], where we use one level of HYP smearing [52 for the temporal gauge 
links in order to improve the signal-to-noise ratio of the correlation functions [53|- For the up, down, and strange 
quarks, we use a domain- wall action [48H50j . and the gluons are implemented using the Iwasaki action |54l 155] . Our 
calculations are based on gauge field configurations generated by the RBC/UKQCD collaboration [51] using these 
actions. Further details of the lattices will be given in Sec. |IIIB[ 

In order to extract the matrix element ([3]), we need to compute suitable three-point- and two-point-functions as 
discussed in the following. We use the following baryon interpolating fields. 
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where a,b,c are color indices, a,/?, 7 are spinor indices, and C is the charge conjugation matrix. The tilde on the 
up, down, and strange quark fields u, d, s denotes gauge-covariant three-dimensional Gaussian smearing, intended to 
reduce excited-state contamination in the correlation functions. 

In the three-point functions, we use the following C'(a)-improved discretization of the continuum HQET current 
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The current J, 



(HQET) 



is renormalized in the MS scheme at = m^. Note that here we match from lattice HQET 



to continuum HQET, but not yet to QCD. This is important because the form factors -Fi and F2 are defined in 
continuum HQET, not full QCD. The matching to QCD will lead to radiative corrections to the simple relationship 
([3]) which depend on F. We will return to this issue in Sec. IV when computing the differential branching fraction for 

In Eq. (11), symmetries of the lattice actions and the equations of motion have been used to reduce the number 



of operators that appear |56j . The term with coefficient Cp provides ©(m^a) improvement, while the term with 

coefficient Cp^""'* provides 0{psa) improvement (here Ps denotes the momentum of the strange quark). The quantity 
Wq^^ is related to tadpole improvement, and is defined as — 1 — aTO5 + 4(l — uq) [SS], where 0777.5 is the domain- wall 

height and Uq is the 4th root of the average plaquette. The matching coefficients Z, c- 



(rriga) 



and 



have been 



computed to one-loop order in lattice perturbation theory for the actions used here in Ref. 56 and are evaluated 
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FIG. 1. Propagator contractions for the forward three-point functions (left) and backward three-point functions (right). The 
thick vertical line at the spatial point y indicates the static heavy-quark propagator. The source for all light and strange-quark 
propagators is located at the fixed point (a;o,x). We sum over all points y, with the appropriate momentum phases as in 
Eqs. (flSl and |T9|. 



at the scale ~ a^^ (the inverse lattice spacing). The coefficient Z is independent of F, but Cp"'"' and c^^""^ 
change sign depending on whether F commutes or anticonimutes with 7*^. We are interested in the matrix element ([3| 
renormaUzed at /x = mf,, and following Ref. (56) we therefore perform a renormalization-group (RG) evolution from 
/i = a^^ to /X = mb, using the two-loop anomalous dimension of the heavy-light current in HQET, which was derived 



in Refs. [57l[58]. This leads to the multiplicative factor U{mb,a~^) in Eq. (11). The RG running is performed with 
Nf = 3 flavors from ^ — a^^ down to ^ — rric, and then with Nf — 4 flavors from /i — rric up to ^ — mi,. This 
two-step running is used because the nonperturbative lattice calculations are done with Nf = 2 + 1 dynamical flavors, 
and with > nic. However, note that doing a simple Nf — 4 running from // = to fi = mi, gives a result that 

differs only by 0.5%. Numerical values for U(mi,,a~^), Z, c^^^°'\ c^'°'\ and ug will be given in Table in the next 
section. 

Having defined the interpolating fields and the current, we will now discuss the correlation functions. We compute 
"forward" and "backward" two-point functions for the A and Aq as follows: 

Ci:^\p',t) = ^e-'P'-(''--) (A5(xo + t,y) A„(a;o,x)) , (14) 
y 

y 

= (Aq,-(xo + <,x) Aq„(xo,x)), (16) 
^(2,A«,bw)^^^ = (Aq,(xo,x) AQ„(xo-t,x)), (17) 



where the superscript "bw" denotes the backward correlator. In Eqs. (16) and (17), the Aq interpolating fields at 
source and sink are required to be at the same spatial point x because of the static heavy-quark propagator. Finally, 
the forward and backward three-point functions for a given gamma matrix F in the current are defined as 

C(;^)(F, p',t,i')-E^"'''''^""''^(^*(^0'^) 4'"^''^^\^o-t + t\y) AQ„(xo-t,y)), (18) 
y 

clf'"\T,p\t,t-t')^J2^''''''^'''^H^Q-i^o + ^^y^ 4«'5^^)(xo + i',y) A5(xo,x)). (19) 



All of the above correlation functions ( 14)-(|19| can be computed using light and strange quark propagators with a 
Gaussian-smeared source located at (a;o,x). For the three-point functions, the quark propagator contractions are 
illustrated schematically in Fig. [l] Because no additional domain-wall propagators are required, we can efficiently 
compute the three-point functions for arbitrary values of t and t', only limited by statistical precision. 

In order to discuss the spectral decomposition of the correlation functions, we introduce the following definitions 
for the overlap factors: 

(0|Aq„(0)|Aq(s)) =Za^Z^„(s), (20) 
(0|A„(0)|A(p',s)) = [(Zf +Zf 70)«(p',s)]„, (21) 
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where Zj^"^ and Z^^-* depend on p'. Here we need two different overlap factors Zj^^ and Zj^"^ for the A, because the 
spatial-only smearing of the quarks in the interpolating field (10 1 breaks hypercubic symmetry |59| . The spectral 



decompositions of the two-point and three-point functions then read 
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where we have only shown the ground-state contributions, and the ellipsis denote excited-state contributions that 
decay exponentially faster with the Euclidean time separations. For the three-point functions, we have used Eq. ([s]) 
to express the current matrix element in terms of the form factors Fi and F2. 
Using the three-point and two-point functions, we then define the following ratio, 



Tl{T,p',t,t') = 



Tr[C(2.A.av)(p,^i)] Tj.[(^(2,AQ,av)(^)] 



(26) 



where the traces are over spinor indices, and the two-point functions in the denominator are the averages of the 
forward- and backward two-point functions (to increase statistics). For the ground-state contributions, the product of 
forward and backward three-point functions in the numerator of Eq. (26) eliminates the i'-dependence, and dividing 



by the two-point functions evaluated at the same t also cancels the t-dependence and Z-factors. For gamma matrices 
r that commute (anticommute) with 7", the ground-state contribution in the ratio TZ{T,p' ,t,t') will be proportional 
to [F+]2 ([F_]2). Thus, we form the combinations 



7^+(p', t, t') = J [7^(l, p', t, t') + n{j^-f^ p', t, t') + 7^(7^^ p', t, t') + nij'j^,p', t, t')] , 
7^_(p',^,^') = ^ [7^(7^p',^,^') + 7^(7^p^^,^') +7^(7^p^^,^') + 7^(75,P^^,^')] , 



which are equal to 



7^+(p',^,^') 



Ea - niA 



n_ip',t,t')^i^^^^[F^]' + ..., 



(27) 
(28) 

(29) 
(30) 



where, as before, the ellipsis denote excited-state contributions. Note that multiplying the gamma matrices in Eqs. ( 27 ) 



and ( 28 ) with 7° would not give any new information, because 7*^(3 = Q. Next, we average ( 27 ) and ( 28 ) over momenta 
p' with fixed magnitude |p'|, and replace the label p' by jp'p to denote the direction-averaged quantities, 



7^±(|p'|^^,0■ 



(31) 



Finally, we evaluate TZ±{\p'\'^ ,t,t') at t' = t/2 (or average it over {t — a)/2 and {t + a)/2 for odd values of t/a) 
where the excited-state contamination is smallest, rescale using i<^A(|p'P) and niA obtained from fits to the two-point 
functions, and take the square root to obtain 
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For t 00, the quantities i?±(|p'| ,t) become equal to the form factors F±{Ea) where Ea = £^a(|p' 
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Set 


N! X Ntx N5 


arris 


(sea) 

arris 


(sea) 


a (fm) 


(val) 

ami 


(val) 


mi^^^ (MeV) 


m^7^ (MeV) 


mcas 


C14 


24^ X 64 X 16 


1.8 


0.04 


0.005 


0.1119(17) 


0.04 


0.001 


245(4) 


761(12) 


2705 


C24 


24^ X 64 X 16 


1.8 


0.04 


0.005 


0.1119(17) 


0.04 


0.002 


270(4) 


761(12) 


2683 


C54 


24^ X 64 X 16 


1.8 


0.04 


0.005 


0.1119(17) 


0.04 


0.005 


336(5) 


761(12) 


2780 


C53 


24^ X 64 X 16 


1.8 


0.04 


0.005 


0.1119(17) 


0.03 


0.005 


336(5) 


665(10) 


1192 


F23 


32^ X 64 X 16 


1.8 


0.03 


0.004 


0.0849(12) 


0.03 


0.002 


227(3) 


747(10) 


1918 


F43 


32-^ X 64 X 16 


1.8 


0.03 


0.004 


0.0849(12) 


0.03 


0.004 


295(4) 


747(10) 


1919 


F63 


32^ X 64 X 16 


1.8 


0.03 


0.006 


0.0848(17) 


0.03 


0.006 


352(7) 


749(14) 


2785 



TABLE I. Parameters of the gauge field ensembles and quark propagators. Here, TVs is the extent of the 5th dimension of the 
lattice, and ams is the domain-wall height [51]. The sea quark masses am,''"^' were used in the generation of the ensembles, and 
we use the valence-quark masses am,™'' when computing domain-wall propagators. The values for the lattice spacings, a, are 
taken from Ref. 16O4 . We denote the valence pion masses by mi^^' , and m^^^' is defined as the mass of the pseudoscalar meson 
with valence strange-antistrange quarks, but without any disconnected contributions (we use mj,^^' to tune the strange-quark 
mass, using the approach of Ref. [61]). Finally, A^moas is the number of light/strange domain- wall propagator pairs computed 
on each ensemble. 



a (fm) 


U{rnb,a ^) 


Ul) 


Z 




„(Psa) 


0.112 


1.09964 


0.875789 


0.9383 


-0.1660 Gr 


-0.1374 Gr 


0.085 


1.06213 


0.885778 


0.9526 


-0.1482 Gr 


-0.1294 Gr 



TABLE II. Renormalization parameters for the matching of LHQET to HQET in the MS scheme, from Ref. [56]. Here, Gr is 
defined by 7"r7" — GrL, so that Gr = +1 if F commutes with 7", and Gr = — 1 if F anticommutes with 7". 



B. Lattice parameters 



The details of the domain- wall/Iwasaki gauge field ensembles generated by the RBC/UKQCD collaboration can 
be found in Ref. [SI]- In Table |l) we summarize the main properties of the subset of ensembles used here, as well as 
some parameters of the domain-wall propagators that we computed on them. There are ensembles with two different 
lattice spacings a « 0.11 fm and a ~ 0.085 fm, with lattice dimensions of 2A^ x 64 and 32'^ x 64, respectively, so that 
the spatial box size is L w 2.7 fm in both cases. We will refer to these two lattice spacings as "coarse" and "fine". 
At the coarse lattice spacing, we use only one ensemble with the lightest available up/down sea-quark masses. At the 
fine lattice spacing, we use two different ensembles. 

In order to construct the correlation functions discussed in Sec. [Ill A[ we require domain-wall propagators with 
Gaussian-smeared sources at (a;o,x), and with masses corresponding to the strange quark as well as the (degenerate) 
up/down quarks. As shown in Table |lj we have seven different combinations of parameters, which we denote as C14, 
C24, C54, C53, F23, F43, F63 (where C, F stand for "coarse" and "fine", and the two digits indicate the hght and 
strange valence quark masses). In four of these combinations, the valence-quark masses are chosen to be lighter than 
the sea-quark masses ( "partially quenched" ) , while the other three combinations have valence-quark masses equal to 
the sea-quark masses (unitary case). On each gauge configuration, we use 0(10) source locations (xcx) to increase 
statistics. The resulting total numbers of "measurements", A^mcasj a-re listed in Table [l] On each configuration, we 
average the correlators over the source locations prior to further analysis. 

In the static heavy-quark action, we use gauge links with one level of HYP smearing with the parameters 
(01,02,0:3) = (1.0, 1.0, 0.5) as introduced in Ref. [62] • The numerical values of the matching coefficients needed for 
the current (111 are taken from Ref. [56 and are given in Table [n] for the choice of HYP smearing parameters used 
here. 



Results for 71+ and TZ- 



At the coarse lattice spacing, we computed the three-point functions (18), (19) for the source-sink separations 
t/a = 4, 5, 15, and at the fine lattice spacing for t/a = 5, 6, 20. We computed these three-point functions 
for lattice momenta p' with < |p'p < 9 • (27r)^/L^. We then constructed the quantities (27) and (28) using 



statistical bootstrap with 1000 samples. When performing the momentum direction average for the largest momentum 
Ip'P = 9-(27r)2/L2^ we used only p' = (2,2, l)-27r/L and lattice symmetries applied to that (for |p'p < 9-(27r)2/L2^ all 
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possible p' with the same magnitude are related by cubic rotations and reflecti ons, a nd we average over all of them). 
Examples of numerical results for the quantities Tl±{\p'\'^ ,t,t') defined in Sec. Ill A are shown in Figs. 2][5 Except 



in the immediate neighborhood of <' = and t' = t, the results for Ti-±{\p' \,t,t'} show only a weak dependence on 
the current-insertion time t' . However, a significant dependence on the source-sink separation t is seen, in particular 
for 7?._. Consequently, we need to extrapolate the results to infinite source-sink separation in o r der to remove the 
excited-state contamination. We perform these extrapolations for i?±(|p'P,t) [defined in Eqs. (32), (33)], as discussed 
in Sec. HITeI 
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FIG. 2. Numerical results for K±{\p'\'^ ,t,t') at |p'p = 1 ■ {2Tv/Lf from the C54 data set (a = 0.112 fm, am 
am^'^^ = 0.005). The source-sink separations shown here are (from left to right) t/a = 6, 8, 10. 
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FIG. 3. Like Fig. [2] but at Ip'l^ = 4 • (27:/ Lf 
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FIG. 4. Numerical results for TZ±{\p' l"^ ,t,t') at |p'p = 1 ■ i2n/Lf from the F43 data set (a = 0.085 fm, am 
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am^'^^ ~ 0.004). The source-sink separations shown here are (from left to right) t/a = 8, 11, 14. 
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FIG. 5. Like Fig.|4] but at \p'f = 4- {2TT/Lf 



D. Fitting the A two-point functions 



In order to compute the quantities R±{\p'\'^ ,t), which were defined in Eqs. (32), (331, we use the energy E'Adp'P) 
and mass ttt-a = E\{Q) of the A baryon computed on the lattice for the same data set. To this end, we average the 
two-point function Tt[C^^-^'^''^ {p' ,t)] over momenta p' with the same |p'P, and perform correlated fits of the form 
A e~^'^^ for sufficiently large t so that the excited-state contamination in the fitted Ea is negligible compared to the 
statistical uncertainty. To give an idea of the signal quality, we show the effective energy (ai?A)cff — \n[C (t) / C {t + a)] 
for selected momenta and two data sets in Fig. [6] The complete fit results for gEa are listed in Table When 
computing R±{\p'\'^,t) via Eqs. (32), (33), we used bootstrap to fully take into account the correlations between E\, 
TUA, and 7^±(|pf , t, t/2). 
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FIG. 6. Effective-energy plots for the A two-point functions at selected momenta. The shaded bands indicate the energies 
extracted from exponential fits of the two-point functions, as well as the t-range used for these fits. Left panel: a — 0.112 fm, 
ami™'' = 0.04, am^jf = 0.005; right panel: a = 0.085 fm, ami™'' = 0.03, am^™'' = 0.004. 
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|p'|V(2^/L)2 


C14 


C24 


C54 


C53 


F23 


F43 


F63 





0.7046(41) 


0.7135(45) 


0.7350(34) 


0.7150(42) 


0.5191(41) 


0.5353(29) 


0.5540(17) 


1 


0.7542(46) 


0.7613(53) 


0.7821(40) 


0.7653(49) 


0.5575(46) 


0.5717(32) 


0.5894(19) 


2 


0.8027(58) 


0.8061(69) 


0.8274(52) 


0.8151(62) 


0.5931(53) 


0.6062(38) 


0.6240(24) 


3 


0.8460(76) 


0.8495(97) 


0.8709(69) 


0.8612(84) 


0.6288(67) 


0.6422(50) 


0.6568(31) 


4 


0.889(10) 


0.896(11) 


0.9106(87) 


0.892(11) 


0.6652(89) 


0.6727(66) 


0.6892(43) 


5 


0.921(11) 


0.929(13) 


0.948(11) 


0.928(13) 


0.695(10) 


0.7033(79) 


0.7201(52) 


6 


0.954(16) 


0.965(19) 


0.989(14) 


0.966(17) 


0.730(13) 


0.739(10) 


0.7479(64) 


8 


1.006(27) 


1.021(31) 


1.051(24) 


1.014(30) 


0.800(23) 


0.802(17) 


0.804(11) 


9 


1.033(35) 


1.066(45) 


1.087(32) 


1.058(40) 


0.839(30) 


0.837(21) 


0.822(14) 



TABLE III. Fit results for aE'Adp'P) from the different data sets (see Table |l]). 



E. Extrapolation of R+ and R- to infinite source-sink separation 



Using the methods outhned in the previous sections, we have obtained numerical results for i?±(|p'P, t) for multiple 
momenta |p'P, for a wide range of source-sink separations t, and for seven different combinations of quark masses 
and lattice spacings. The next step is to compute the ground-state form factors F± by extrapolating R± to infinite 
source-sink separation. In the following, we use the notation R''^{t), where i = C14,C24, ...,F63 labels the data set 
(see Table|l|, and n = 0, 1, 2, 3, 4, 5, 6, 8, 9 labels the momentum of the A, writing |p'p = n ■ (27r)^/i^. 

At zero momentum, we can only compute i?+ because TZ^ vanishes for = m/^ [see Eq. mA]. Results for i?!^"(i) 
from the two data sets i = C54, F43 are plotted in Fig. [7] as a function of the source-sink separation t. Remarkably, 
no significant t-dependence is found in this quantity, allowing constant fits of the form 



i?f(t) 



(34) 



which are also shown in Fig. [t] The fits fully take into account correlations and have x^/dof < 1, showing that there 
is indeed no evidence for deviation from a constant. Note that i?+ can be obtained from the vector current sj'^b (in 
our calculation we have replaced 7" by 1 because j'^Q = Q for static heavy quarks). In relativistic QCD, neglecting 
mass effects, charge conservation would then prevent any t-dependence at zero momentum, because excited states 
have the same charge as ground states. It appears that some remnant of this symmetry also remains in our case. 

At non-zero A-momentum, both R^ and i?_ can be extracted, and significant t-dependence is seen. To extract the 
ground-state contributions and F-, we perform fits including an exponential correction term that describes the 





t (fm) 



FIG. 7. Constant fits to the t-dependence of R+{t) at |p'p — 0. Left panel: 



0.112 fm, am. 



(val) 



0.04, am^™'' = 0.005; 



right panel: a = 0.085 fm, ami™''' ~ 0.03, am}^^^' = 0.004. The fits have x'^/dof values of 0.92 and 0.78, respectively. 



(val) 



10 



leading effects of excited states, 



R±. [t) ^ exp[-d_(. t\, 



(35) 



where _F!|.'", AljJ^, and i5^" are the fit parameters, which exphcitly depend on the data set i and the momentum n. 
Because the energy gaps are positive by definition, we write 



Sr/il GeV) = exp(Zr), 



(36) 



and use instead of i5!|!" as the fit parameters. The fits using Eq. (36 1 are performed separately for each momentum 



n, but simultaneously for the different data sets i. Note that the size of the momentum unit, {2tt)/ L (in GeV), is the 
same at the coarse and fine lattice spacings within uncertainties, because the box sizes (in physical units) are equal 
within uncertainties. Performing the fits simultaneously for the different data sets at the same momentum allows us 
to use the prior knowledge that the hadron spectrum does not change dramatically when the lattice spacing or quark 
masses are varied by small amounts. To this end, we augment the function used to perform the fits to Eq. (35 1 as 
follows: 



/iC14,n 

2 I 

X +- 



C14,C24-i2 
m J 

. _ ;F43,„^ 



F23,F43i2 



/7C24.n 



C54,nN,2 



C24X5412 



/7F43.ri 



F63,n\2 
± ) 



F43,F63i2 



C53,nN2 
± } 



r C54,C53i9 
/iC54.n 7F63.nN2 



C54,F63i2 



(37) 



where 



i \2 



(mi)1^+<[(m;j^-(m; 



\2l2 



(38) 



with Wra = 4 GeV , and (Tq — 0.1. With these parameters, Eq. (371 implements the constraint that the energy gaps. 



at given A-momentum n, should not change by more than 10% when going from the fine to the coarse lattice spacing, 
and not more than 400% times the change in or ml^ (in GeV^). Note that absolute variations of l^J^ translate to 



relative variations of , because d[exp(l 



2,n \ 
± 



dZ^". 



Example fits of (t) using Eq. ( 35 ) are shown in Fig. [sj The fits are fully correlated, using covariance matrices 
computed from the bootstrap ensembles for RjJ^{t). Note that the excited-state contribution in which is negligible 
at p' = (cf. Fig. [t]), gradually increases with the momentum. In contrast, i?„ shows the strongest excited-state 
overlap at the smallest momentum, and this overlap decreases as the momentum increases. The excited-state overlap 
is slightly stronger at the fine lattice spacing when compared to the coarse lattice spacing. This is expected because 
the quark smearing width in the baryon operators was different for the two lattice spacings (we used the same width 
in lattice units). We only computed the correlators for </a > 4 at the coarse lattice spacing and t/a > 5 at the fine 
lattice spacing. At the fine lattice spacing, it was necessary to exclude the points with t/a < 8 from the fits to R-. 
Once these points were excluded, all fits had x^/dof ss 1.0. Given the limited time range and the limited statistical 
precision of the available data, it was not possible (and not necessary) to perform fits with more than one exponential. 
As a check, we have also performed fits without the constraints (37 1, which give consistent results but are less stable. 

The fitted values of the energy gap parameters, (5^" — exp(Z!|l")-(lGeV), are shown as a function of the A-momentum 
for one ensemble in Fig. [9] (left panel). Within uncertainties, we find that 



0+ 



(39) 



for all data sets i and momenta n. The energy sp ectr um is a property of the QCD Hamiltonian and is independent of 
the correlation function considered, so the result (39 1 is not surprising. However, one possible situation in which 
and S'l^ would be different is when an excited state has negligible overlap in but significant overlap in i?_ (or vice 
versa). Furthermore, by using only a single exponential, we may be effectively averaging over multiple excited states 
which we cannot resolve individually, but which may appear with different sets of weights in and i?_ . Having said 
that, the values of Sl^" and 5!^" from our fits are in complete agreement and it is evident that the separate parameters 

and (5!:" may be replaced by a single parameter (5''". Thus, we performed new, coupled fits of R+ and i?_ of the 
form 



R^"{t) = Fi" + exp[-,5''" t], 



(40) 
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FIG. 8. Fits of the i-dependence of R±{t) using Eq. (|35|. Left panels: a = 0.112 fm, am'i""'' = 0.04, am^™'^ = 0.005; right 



3.0 



panels: a = 0.085 fm, ami™'' ~ 0.03, am^™'' = 0.004. At the fine lattice spacing, only the points with t > 0.6 fm are included 
in the fit of R-. The maximum values of t for the data were limited by statistical noise in the two-point and three-point 
functions; for t larger than the right-most points in the plots, the statistical fluctuations were too large to compute the square 
roots in Eqs. ((32|| and pSl). 
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FIG. 9. Left: results for the energy gap parameters from separate fits of and 7?_ using Eq. (351. The points are offset 
horizontally for clarity. Right: results for the energy gap parameters from coupled fits of R+ and R- using Eq. (401. The data 

-0.04, am^™" = 0.005. 



shown here are for a = 0.112 fm, arris 



|pf/(2^/L)2 


C14 


C24 


C54 


C53 


F23 


F43 


F63 





0.9416(56) 


0.9443(66) 


0.9430(59) 


0.9390(71) 


0.9320(58) 


0.9281(49) 


0.920(13) 


1 


0.868(14) 


0.874(27) 


0.871(23) 


0.873(26) 


0.878(33) 


0.876(18) 


0.807(25) 


2 


0.785(15) 


0.786(34) 


0.788(27) 


0.791(25) 


0.823(25) 


0.818(14) 


0.738(28) 


3 


0.717(15) 


0.705(36) 


0.715(28) 


0.723(24) 


0.770(26) 


0.763(14) 


0.680(28) 


4 


0.671(15) 


0.669(26) 


0.672(22) 


0.668(17) 


0.697(32) 


0.698(15) 


0.633(28) 


5 


0.618(15) 


0.612(26) 


0.621(21) 


0.622(15) 


0.667(27) 


0.664(16) 


0.587(27) 


6 


0.586(16) 


0.580(20) 


0.576(20) 


0.577(15) 


0.625(30) 


0.621(19) 


0.554(29) 


8 


0.537(14) 


0.545(16) 


0.532(13) 


0.516(24) 


0.569(27) 


0.563(22) 


0.507(39) 


9 


0.508(17) 


0.516(16) 


0.502(15) 


0.491(26) 


0.546(31) 


0.536(28) 


0.484(40) 



TABLE IV. Results for _F+ from the different data sets (see Table [T|. The uncertainties presented combine statistical and 
systematic fitting uncertainties in quadrature. 



with energy gap parameters 5*'" — exp(P'") • (1 GeV) that are shared between and The results for 6^'" from 
the coupled fits are shown in the right panel of Fig. m We note that the coupled fits had values of x^/dof that were 
as good or better than the values from the separate fits, confirming that the assumption (39 1 is justified. 

The results for the extracted ground-state form factors F^'" from separate and coupled fits are compared in Fig. 10 
As can be seen there, the results are consistent with each other, but the coupled fits with shared energy gap parameters 
give significantly smaller uncertainties, in particular for Therefore, we use the results from the coupled fits in our 
further analysis. 

To estimate the systematic uncertainty resulting from the choice of the fit range in t, we computed the changes in 
F^l"' when increasing t^i^ by one unit simultaneously for all data sets, thereby removing the points with the most 



contamination from additional higher excited states. As can be seen in Fig. [TT] the resulting change in the fitted 
F^"' is small, and in most cases consistent with zero. Nevertheless, we add this shift in quadrature to the original 
statistical uncertainty of F^" . The final results for F^" , including this systematic uncertainty, are given in Tables 
and[v} For convenience, we also provide results for Fl'^ — {F^" + F'2^)/2 and Fj'" = (F|'" — Fl'")/2, computed 



IV 



using bootstrap to take into account the correlations between F^ and Fj , in Tables 



VI 



and 



VII 
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FIG. 10. Fit results for the ground-state contributions _F+ and from separate fits of i?+ and _R_ using Eq. (|35l) vs results 



arris 



■ 0.03, am^' 



(val) 



from coupled fits using Eq. (40 1. Left panel: a — 0.112 fm, ami™'' ~ 0.04, am'i^'^^' = 0.005; right panel: a — 0.085 fm, 
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0.004. The points are offset horizontally for clarity. 
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FIG. 11. Fit results for the ground-state contributions F+ and F- from coupled fits of 7?+ and R- using Eq. (401: effect of 



increasing all tmin values by one unit. Left panel: a = 0.112 fm, ami™'' = 0.04, am^™' = 0.005; right panel: a — 0.085 fm 
ami™'' = 0.03, am^J^'^ = 0.004. The points are offset horizontally for clarity. 
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C14 


C24 


C54 


C53 


F23 


F43 


F63 


1 


1.355(53) 


1.422(52) 


1.379(33) 


1.346(41) 


1.352(57) 


1.339(60) 


1.296(51) 


2 


1.233(58) 


1.278(36) 


1.274(26) 


1.230(35) 


1.241(45) 


1.258(38) 


1.185(49) 


3 


1.143(36) 


1.163(43) 


1.159(41) 


1.135(72) 


1.183(52) 


1.172(42) 


1.094(47) 


4 


1.032(51) 


1.084(40) 


1.089(29) 


1.071(44) 


1.094(81) 


1.088(52) 


1.004(49) 


5 


0.973(34) 


0.990(35) 


0.990(28) 


0.967(49) 


1.035(73) 


1.029(45) 


0.928(47) 


6 


0.904(35) 


0.905(43) 


0.922(42) 


0.907(65) 


0.938(82) 


0.940(44) 


0.872(45) 


8 


0.752(33) 


0.790(42) 


0.801(37) 


0.820(74) 


0.805(58) 


0.826(33) 


0.742(48) 


9 


0.700(43) 


0.706(51) 


0.719(52) 


0.712(76) 


0.748(58) 


0.778(37) 


0.707(41) 



TABLE V. Results for _F_ from the different data sets (see Table |T|. The uncertainties presented combine statistical and 
systematic fitting uncertainties in quadrature. 



|pf/(2Vi)' 


C14 


C24 


C54 


C53 


F23 


F43 


F63 


1 


1.111(31) 


1.148(32) 


1.125(18) 


1.110(23) 


1.115(40) 


1.107(38) 


1.051(35) 


2 


1.009(35) 


1.032(26) 


1.031(23) 


1.011(26) 


1.032(30) 


1.038(24) 


0.961(37) 


3 


0.930(21) 


0.934(35) 


0.937(32) 


0.929(46) 


0.976(33) 


0.968(25) 


0.887(36) 


4 


0.852(29) 


0.876(20) 


0.881(17) 


0.870(25) 


0.895(50) 


0.893(28) 


0.818(35) 


5 


0.795(18) 


0.801(22) 


0.805(18) 


0.795(27) 


0.851(41) 


0.846(22) 


0.757(34) 


6 


0.745(19) 


0.743(26) 


0.749(26) 


0.742(33) 


0.782(48) 


0.780(24) 


0.713(34) 


8 


0.644(17) 


0.667(24) 


0.666(19) 


0.668(31) 


0.687(32) 


0.695(17) 


0.624(41) 


9 


0.604(23) 


0.611(27) 


0.611(24) 


0.602(32) 


0.647(35) 


0.657(19) 


0.595(35) 



TABLE VL Results for Fi from the different data sets (see Table The uncertainties presented combine statistical and 
systematic fitting uncertainties in quadrature. 



Ipf /(2^/L)2 


C14 


C24 


C54 


C53 


F23 


F43 


F63 


1 


-0.244(23) 


-0.274(26) 


-0.254(23) 


-0.236(25) 


-0.237(23) 


-0.232(23) 


-0.244(19) 


2 


-0.224(25) 


-0.246(23) 


-0.243(14) 


-0.220(16) 


-0.209(21) 


-0.220(16) 


-0.224(16) 


3 


-0.213(17) 


-0.229(20) 


-0.222(13) 


-0.206(27) 


-0.207(24) 


-0.205(19) 


-0.207(15) 


4 


-0.181(24) 


-0.208(27) 


-0.209(19) 


-0.202(22) 


-0.199(35) 


-0.195(26) 


-0.186(18) 


5 


-0.177(19) 


-0.189(21) 


-0.184(17) 


-0.173(24) 


-0.184(37) 


-0.183(26) 


-0.171(18) 


6 


-0.159(20) 


-0.163(21) 


-0.173(20) 


-0.165(34) 


-0.157(38) 


-0.159(24) 


-0.159(17) 


8 


-0.108(19) 


-0.123(21) 


-0.135(21) 


-0.152(46) 


-0.118(32) 


-0.131(22) 


-0.118(16) 


9 


-0.096(23) 


-0.095(27) 


-0.109(29) 


-0.111(47) 


-0.101(30) 


-0.121(27) 


-0.111(20) 



TABLE VIL Results for F2 from the different data sets (see Table |T|. The uncertainties presented combine statistical and 
systematic fitting uncertainties in quadrature. 



15 



F. Chiral and continuum extrapolation of the form factors 



The last step in our analysis of the lattice data is to fit the dependence of F^'" on the quark masses, the lattice 
spacing, and on i?A- The form of the dependence is unknown; low-energy effective field theory combining heavy-baryon 
chiral perturbation theory for the A sector and heavy-hadron chiral perturbation theory for the Ag sector may be 
useful over some range of i?A, but not in the region with |p'| > A;^, where A^ ~ 1 GeV is the chiral symmetry-breaking 
scale. We therefore use a simple ansatz that fits our data well at the level of statistical precision that we have. In the 
following it is advantageous to express the form factors as functions of the energy difference Ej^ — mj^ instead of i?A, 
as this depends less on the quark masses. We find that this dependence can be described well using a dipole function 
of the form F± — N±/{X± + E\ — miCf' . We generalize this ansatz to allow for dependence on the light and strange 
quark masses, as well as the lattice spacing, in the following way: 



iV4 



{XI 



where the functions X\. are defined as 



X\ 



X± + Q,± • [{mir - {m^'^n'] + Cs,± ■ 



(41) 



(42) 



As before, we use the notation where i = C14, C24, F63 labels the data set (see Table|l]), and n labels the momentum 
of the A. The free fit parameters in Eq. (41 1 are iV±, X±^ d±, ci^±, and Cs^±. The dependence of the form factors on 



the light and strange quark masses is described by allowing X^ to depend linearly on {m]^) and (m^_^) , where 



and m^^ are the valence tt and rjs masses for each data set i, as given in Table |l] We wrote Eq. (42) in terms of the 
differences between the lattice and physical masses for convenience, with m^^^^ — 138 MeV and mf^^^ = 686 MeV 
|61j . The leading dependence of the form factors on the lattice spacing is expected to be quadratic in a, owing to the 
chiral symmetry of the domain- wall action and the use of the order-a- improved current ( 11 ). Discretization errors are 



expected to grow as the mom entum of the A increases. We therefore incorporate the a-dependence using the factor 



l + d±{a'EY"n in Eq. (41) 



In our fits, we take into account the correlations between the results for F^"' at different momenta 



n and different 

data sets i (in the case where the data sets correspond to the same underlying ensemble of gauge fields). The fits 
are performed independently for F^"' and F^" . To account for the uncertainties and correlations of the A baryon 
energies E^ 

-ii 



and I' J 

(including the masses m\ = E^j^) in Eq. (41), we promote E]'"' 



to additional parameters of the fit, and 



add the term Y.i,n,i' ,n' [Cov(Ea) ]i,n,i',n' {E\ 



A ' 



Ejy' ) to the X function, where Ejl are the previous 



results from the fits to the two-point functions, and the energy correlation matrix Cov(£'a) was computed from the 
bootstrap ensemble of the two-point fit results. Using a similar term, we investigated the inclusion of the further 
correlations between the A energies and the form factor values F^^\ but with the current level of statistics, such fits 
did not converge to a stable minimum of x^- 
The fits using Eq. 



Fi" from Tables 



and pseudoscalar masses, 



41 ) are visualized as a function of E\ — m\ in Fig. 12 There, we show the results for 



V| along with the fitted functions (41) evaluated at the corresponding lattice spacings a* 
and m]^^ . The data are described well by the fitted functions (the F63 set fluctuates 

shows the fit functions 
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downward, but the overall values of x^/dof are smaller than 1). The bottom-right plot in Fig, 
evaluated in the continuum limit (a — 0) and for the physical values of the pseudoscalar masses. By construction, in 
this physical limit, Eq. (41) reduces to 



- (X,+FA-mA)- 

which only depends on the parameters N± and X± . Our results for N± and X± are given in Table |VIII[ The results 
for the parameters describing the dependence on the quark masses and the lattice-spacing are c;.+ — 0.094(32) GeV~^, 
Cs,+ = -0.019(27) GeV"\ d+ = 0.027(27), q,_ = 0.04(20) GeV"\ c^,- = -0.14(11) GeV"\ and d_ = -0.036(67), 
which are all very small and mostly consistent with zero. 

Functions for the form factors Fi and F2 could be obtained from ( [43| by taking the linear combinations Fi — 
{F+ +F-)/2 and F2 = (F+ — F_)/2. However, because we use independent pole parameters X+ and X_, these linear 
combinations are no longer of the simple dipole form. Alternatively, we can also perform new fits to the lattice data 
pi,n _ ^phn _j_ p^^'^y2 and F2'" = (F^'" — F!:")/2 using functions of the same form as in Eq. (41), but with new 
parameters labeled by the subscripts 1, 2 instead of +, — : 

Fl:2 = V—^n ■ [1 + d,.2ia^Ern (44) 



F 
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Parameter 


Result 


N+ 


3.188 ±0.268 GeV^ 


x+ 


1.852 ±0.074 GeV 


N- 


4.124 ±0.750 GeV^ 




1.634 ±0.144 GeV 



TABLE VIII. Fit results for iV± and X± using Eg. ( |4l[ ). The covariances are Cov{N+,X^ 
0.106 GeV^ The results are renormalized in the MS scheme at ^ = rrib. 



= 0.0198 GeV^ Cov{N-,X-) = 



Parameter 


Result 


Ni 


3.975 ± 0.576 GeV^ 


Xi 


1.776 ±0.123 GeV 


N2 


-0.385 ±0.132 GeV^ 


X2 


1.156 ±0.200 GeV 



TABLE IX. Fit results for A^i,2 and Xi^2 as discussed in the main text. The covariances are Cov{Ni, Xi) — 0.0692 GeV^, 
Cov(A'^2,^2) = —0.0256 GeV^. The results are renormalized in the MS scheme at /j, — mt. 



These fits are visualized in Fig. 13 and the resulting parameters A^i,2, ^1,2 are given in Table IX In this case, 
the results for the other fit parameters were q^i = 0.09(17) GeV"\ c^.i = -0.067(94) GeV~\ di = -0.049(53), 
= -0.06(38) GeV-\ Cs,2 = -0.35(22) GeV~\ and da = 0.00(15). 



G. Estimates of systematic uncertainties 



The remaining systematic uncertainties in our form factor results include missing higher-order renormalization 
corrections to the heavy-light current, finite-volume effects, chiral-extrapolation errors, and residual discretization 
errors. We discuss each of these below. Furthermore, for large Ea — rriA, where we do not have lattice data, our 
assumption of a dipole shape in Eqs. ( 41 ) and ( 44 ) introduces an unknown model-dependence. This is illustrated in 



Fig. [14j), where we compare the dipole fits to monopole fits. However, we do not have confidence that this difference 
is a reliable estimate of a fitting form systematic uncertainty (or indeed that such a systematic uncertainty can be 
constructed) and so leave this to the judgment of the reader. 

To est ima te the systematic uncertainty due to missing higher-order renormalization corrections to the heavy-light 
current (11), we vary the scale fi in the matching coefficients Z{fi), c^™''°-'(/i), c^p=°-'(/^), and in the renormalization- 
group running U{mi,,fi). We then recompute the ratios (32) and ( |33[ ) with the modified current. Changing fi from 

to 2a~^ results in a 7% change of both i?+ and i?_ at the coarse lattice spacing and a 6% change of both 
an i?_ at the fine lattice spacing. These relative changes are nearly independent of the source-sink separation, the 
momentum, and the quark masses. Thus, we take the renormalization uncertainty in the final form factor results to 
be 6%. 

Finite-volume effects in the lattice data are unknown (as in the chiral extrapolation, no low-energy effective theory 
exists to guide us over the full range of ^a), but are expected to be of order exp(— m^L). The lowest pion mass 
used in our calculation is « 227 MeV, corresponding to mTri « 3.1 and exp(— m^L) « 0.04. The average value 
of exp(— TTiTri) for the different data sets (see Table |l| is about 0.02. Given these values, we estimate the systematic 
uncertainty in our final results due to finite- volume effects to be 3%. 

The chiral extrapolations of the form factors were performed quadratically in the valence pseudoscalar masses, 
i.e. linearly in the valence-quark masses, ignoring that some of the data were partially quenched and ignoring possible 
nonanalytic dependence on the quark masses. To study the effect of the quark-mass extrapolations, we perform new 
fits with with either ci ± or Cs^±, or both, set to zero, and consider the changes in the extracted form factors and 
F_ (analogously also for Fi and F2). This corresponds to replacing the linear fits of the quark-mass dependence by 
constant fits. The resulting relative changes in F+, F-, Fi, and F2 when setting c; = are below 1% throughout the 
kinematic range where we have lattice data; the biggest relative change (5%) is seen in F2 at zero recoil when setting 
Cg — 0. However, all of the changes are consistent with zero within statistical uncertainties. 
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FIG. 12. Fits of the form factor data for F+ and F- using Eq. The fit of F+ has xVdof = 0.84, and the fit of F- has 

X^/dof = 0.72. The lowest right plot shows the continuum form factors for the physical light- and strange quark masses. 
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FIG. 13. Fits of the form factor data for Fi and F2 using Eq. (|44|. The fit of Fi has x^/dof = 0.59, and the fit of F2 has 
X'^/dof = 0.85. The lowest right plot shows the continuum form factors for the physical light- and strange quark masses. 
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FIG. 14. Comparison of our results for the Aq — >■ A form factors from dipole fits using Eqs. (41 1 and (44 1 to results from 
monopole fits (the same equations without the power of 2 in the denominator). Shown here is the entire kinematic range needed 
for the decay A;, — >■ M^l~ with mi = (the point = corresponds to — m\ « 1.8 GeV). In the range where we have 
lattice data [Ea — rriA < 0.7 GeV), the dipole and monopole functions are consistent with each other, but for large Ea — "lA, 
model-dependence can be seen. In our main analysis we choose the dipole fits as they have slightly lower values of x^/dof. 




FIG. 15. Final results for the Aq_h> A form factors, given by F± = N±/{X± +Ea~ to-a)^ and Fi,2 = A''i,2/(A'i,2 + Ea- niAf 
with parameters as in Tables [VIII| and |IX[ The dark shaded bands show the statistical uncertainty, and the light shaded bands 
show the total (including 8% systematic) uncertainty. The results are renormalized in the MS scheme at /x = rrib. 
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FIG. 16. Final results for the form factor ratio ~F2/ Fi, given by -~{N2/N\){X\ + Ea — ?tja)^/ {X2 + EA — mA)^ with parameters 
as in Table [IX] The shaded band shows total uncertainty, which is dominated by the statistical uncertainty. 
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By "residual discretization errors" we mean discretizatio n er rors that are not eliminated th rou gh the continuum 
extrapolation using the factors of [1 + d±ia'-EY')^] in Eq. (4l]) and [1 + di,2(a*-E;^")^] in Eq. (|44|. While we know 
that the leading discretization orders are quadratic in a, we do not know how they depend on E\. To study the effect 
of the factors [1 + d± (a'i?)^")^] and [1 + c?i,2(a*£'^")^], we perform new fits with d± (or ^1^2) set to zero. In this case, 
the resulting changes in all form factors are less than 4% for Ea — toa < 0.8 GeV, and are consistent with zero within 
the statistical uncertainties. 

Combining the uncertainties in the above discussion, we estimate the total systematic uncertainty of our final results 
for the form factors for Ea — itia < 0.8 GeV to be 8%. Plots of the form factors including this systematic uncertainty 



are shown in Fig. 15 The relatively large systematic uncertainty from the current matching cancels in ratios such as 
F2/F1. This ratio is shown in Fig. 16 and we estimate the total systematic uncertainty in F2/F1 to be 5%. 



IV. THE DECAY A;, A £+£' 

As a first application of our results for the Ag — >■ A form factors, we calculate the differential branching fraction for 
the decays A?, — >■ A £^£~ with ^ = e, /z, r in the standard model. Long-distance contributions (discussed at the end of 
this section) are not included. In the following, it is convenient to use the notation 



U 



cff 



i=7,9,10 



(45) 



with 



•^9 



167r2 
16^ 



IGtt- 



167r2 



IGtt' 



and 



(1^ - 

^Y,cB = C'i.cff + C- cff , 



167r2 



(46) 



(47) 



The "effective" Wilson coefficients Ci^cS and C^'^ff (i = 7,9,10), which are defined at the scale /i = rrif,, contain the 
one-loop matrix elements of the four-quark operators Oi,...,Oq 
The invariant matrix element of "Hcff is given by 

M = -{K{p\s') £+ ip+ , s+ ) r (p_ , s_ ) iHeff I Ab (p, s)) 

Gf Ofem 1 



with the hadronic matrix elements 



u{p+,s+)Y'v{p-,S-) + u(p+,s+)7^75w(p_,s_) 



2m6 ^_ 



B„ 



(48) 



Here we have contracted the electromagnetic field strength tensor in Oy' and with a perturbative insertion of the 
leptonic QED interaction term. 

In order to use the HQET relation ([s]) to express the hadronic matrix elements in terms of the form factors Fi and 
F2, we first need to match the QCD b s currents in Eq. (48) to HQET currents. We use the one-loop perturbative 
results in naive dimensional regularization from Ref . [63] : 

sjf^b = S7^Q -I- Cy sVf^Q, 

50-^1/75^ = Co- Sdf^ulbQ: (49) 
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with 



Ca = 1 - 



TT 

3 TT ' 



In 



TT 



"4 


I'M 


-)1 


.3^ 







(50) 



This gives 



2TOf, ^. 2mt ^. \ . \ 



C 



V 



ClO,offC7 7m75 + C^O^cSCv Vf,j5^u{p, s). (51) 



Note that here we use spinors with the standard relativistic normahzation for all particles, including the A^. In terms 
of the HQET spinors (|7|), we have u{p, s) = y/fnX^ U{v, s), with p — m\^v. For a given value of , the form factors 
Fi and F2 in Eq. (51 1 are evaluated at 



Ea^p' -v 



(52) 



where the masses take their physical values. The fully differential decay rate with polarized particles is given by 

1 d^p' d^p_ d^p 



dr 



2mA, (27r)32£;A (27r)32£;^- {2TT)^2Ee 



-(2^)V(p-p'-p_-p+)|X|^ 



(53) 



For the standard model calculation, we set the right-handed couplings to zero (C^ = C'^ = C[q = 0) and use 
the following Wilson coefficients (at ^ = 4.8 GeV), which are of next-to-next-to- leading-logarithm accuracy 

C7,eff = -0.304, 

C9,eff('?') = 4.211 + r(q2), 

Cio,cff = -4.103. (54) 

The function Y{q'^) is defined as in RefJM]- Furthermore, we use \Vts\ = 0.04002 and \Vtb\ = 0.999142 from Ref. [55]. 

To calculate dV/dq^, we integrate (53) over the lepton momenta and the direction of the A, sum over the spins of 
the A, £~ , and average over the A^ spin (because dT/dq^ is rotationally symmetric, it has to be independent of 
the Ah polarization, and therefore we can treat the Af, as unpolarized here). The result is given by 



dr _ a',^Gl\VtbVt 



* 12 

ts\ 



'1 



4mf 



V {{mA^ - ttiaY - q'^){{mA^ + niA)^ ~ q^) 



with 

-4io,io 

-^9,9 
A73 



dg2 6144 TT^ g4 ^5 ^ w 

X g'|Cio,cfifAo,io + 16c>g(g2 + 2mf)\Cj^,s?A7j + q\q^ + 2m2)|C9,cff (g')l' ^9 
+8g^c^TOh(g^ + 2mf)mA,^[C7^cS C9,cff (g^)] A,9 , 

(24 + 2c^c, + cl) {2mf + q^) (m\^ - 2ml^ml + {q^ - mi)') 
+2ml^q^ (4c2 (g2 _ Amf) - (2c-,c„ + c^,) {q^ - lOmf)) ]f + Ac^ (c^ + c„) (2m? -|- q^) QF+F^, 
+ (9' - 2mi) + (g2 - mi)') F + 2gF+F_, 



(2c; + 2c^c„ 



+ cl) (m\^ + (g2 - mlf) - 2m\^ {2cl^ {vfi\ - 2?') + (2c^c„ + c^) (m^ + q^)) 



+4c^ (c^ -I- c^)gF+F_, 



3a 



{ml, -ml + q^)F + 2 (3c^ + c,) (ml - 2ml {K, + + - "^aJ') 



(55) 



(56) 
(57) 

(58) 
(59) 



22 




3.0 

> 

O 2.5 



^ 2.0 

^ 1-5 
a. 

\ 1.0 



< 
t 



0.5 
0.0 







— This work 
^^CDF 2012 



10 

(GeV^) 



15 



20 



FIG. 17. Left panel: Differential branching fraction for At, Afi'^/i^ . The solid curve is our prediction using the form factors 
from lattice QCD. Long-distance effects are not included in the calculation. The inner, dark shaded band around the curve 
indicates the uncertainty in dS/dg^ that results from the statistical plus systematic uncertainty in the form factors F±. The 
outer, light shaded band additionally includes an estimate of the systematic uncertainty in dB/dq^ that results from our use of 
the static approximation for the 6 quark. The vertical dashed line indicates the lowest value of where we have lattice data; 
to the left of that line the form factors are extrapolated. To illustrate the model-dependence resulting from the extrapolation 
of the form factors to low g^, the dashed curve shows dB/dq'^ computed with form factors extrapolated using a different ansatz 
(monopole instead of dipole. 



14 



the uncertainty for the dashed curve is not shown for clarity). The experimental 
[10] . The error bars shown for the experimental data include systematic 
The vertical shaded bands indicate the charmonium veto regions, where long-distance effects are large. Right 



see Fig 

data are from Ref. [B7], which is an update of Ref 
uncertainties 

panel: with binning applied to the theory prediction 



where 



T = ((mA„ - rrif^f - q^)Fl + {{m^^ + m^)^ - q^)Fl, 



3toa 



miy 



(60) 
(61) 



To obtain the differential branching fraction dB/dq 



1.425 • 10- 



TA^dr/dg^, we use the experimental value of the lifetime, 
The form factors F4. and f _ are given by the functions ( 43 ) with parameters N± and X± 

The resulting differential 



as in Table VIII and with additional systematic uncertainties of 
branching fraction for A;, — >■ KjX^ jjT is shown in Fig. 
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included (see Fig. 15). 
along with recent experimental results from CDF 



The agreement of the standard model with the experimental data is clear, with no evid ence for physics beyond the 
standard model. 
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Further predictions for Af, — > Kl'^i with £ = e,T are shown in Fig. 

In Figs. [17] and [18] the inner shaded bands around the curves correspond to the statistical plus systematic uncertainty 
in the form factors F±. However, note that we have lattice data only in the region > 13 GeV^, as indicated by 
the vertical dashed lines in Figs. [17] and [18] Below that region, we rely on extrapolations of the form factors, which 
are model-dependent. This was shown in Fig. 14 where we compared the form factors from dipole and monopole 
fits. Our main results for the differential branching fractions are based on the dipole form factors. To illustrate the 
model-dependence, the dashed curves in Figs. [17] and [TS] give the differential branching fractions calculated with the 
monopole form factors (the uncertainties of the dashed curves are not shown for clarity, but are of similar size as 
with the dipole form factors). In the large-g^ region, both curves are consistent with each other. At low q^, model- 
dependence can be seen, but as already discussed in Sec. |III G[ a comparison between any two fit models can only 
give a qualitative picture of the model-dependence. 

include an estimate of the systematic uncertainty in dB/dq^ which 



The outer shaded bands in Figs. 17 and 
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arises from the use of the static approximation (i.e., leading-order HQET) for the b quark. In general, the uncertainty 
associated with this approximation is of order AqcD/m^. However, the non-zero momentum p' of the A baryon in 
the Afc rest frame is an additional relevant scale, which may lead to errors of order |p'|/mf,. Thus, we add these two 
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FIG. 18. Differential branching fractions for At 
explanations. 



Ae'^e (left) and At — >■ Ar^r (right). See the caption of Fig. 
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for 



errors in quadrature and estimate the relative systematic uncertainty in dB/dq'^ due to the use of HQET to be 



(62) 



where we take Aqcd — 500 MeV. 

Another major cause of systematic uncertainty in the branching fraction is that we have neglected long-distance 
effects. The most important type of long-distance effects is associated with photon exchange between the lepton and 
quark electromagnetic currents, where the quark electromagnetic current combines with the four-quark operators Oi 
to Oe in the effective Hamiltonian ([T]). The resulting contribution to the decay amplitude is described by nonlocal 
hadronic matrix elements of the form 



9 



d^x e"^-(A(p',s')|TO.(0) j,V(:r) |A,(p,s)), (63) 



where j^^^ {x) is the quark electromagnetic current (see for example |681I69| ). These are challenging matrix elements to 
compute from lattice QCD, but may give large contributions when coincides with the mass of a hadronic resonance 
with J^^ = 1~. Of the four-quark operators, Oi = ic''-/''PLb''){s''-ff,PLc'') and O2 = {c''-^''PLb''){s''-/f,PLc'') (where 
the superscripts a and b are color indices) have the largest Wilson coefficients, and at = Wjy^,m^, the decay 
amplitude is dominated by the non-local matrix elements of Oi and 02- The experimental analysis of Ref. [ST] 



excludes these q^ regions, as shown by the vertical shaded bands in Fig. 17 Away from these resonances, however, 
we expect the short distance contributions from 07,9,10 to be the most important. Several approaches show that for 
large g^, these and other long-distance effects can be treated as small corrections to the leading order behavior given 
by matrix elements of 07,9,10 [MMZO]- Although these papers discuss B — ^ K^*'>£^£^, the same principles apply to 



V. CONCLUSIONS 



Theoretical studies of the rare baryon decays Aft — > Ai'^i" and Af, — )■ A7 require knowledge of the hadronic 
matrix elements (A| sTb |Af,) in nonperturbative QCD. At leading order in heavy-quark effective theory, these matrix 
elements are given by two independent form factors, Fi and F2 (or, equivalently, = Fi + F2 and F_ = Fi — F2), 
which are functions of the energy of the A baryon in the A;, rest frame |33H35j . Here, we have performed the first 
lattice QCD calculation of these form factors. Our final results for F± and ^1,2, in the continuum limit and for 
the physical values of the up, down, and strange-quark masses, are shown in Fig. |15| High prec ision determinations 
were achieved by analyzing the ratios _R+(|p'p,t) and i?_(|p'p,t), defined in Eqs. (|32[) and (331, for a wide range of 
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source-sink separations t, and by using multiple light-quark masses as well as two different lattice spacings. Systematic 
uncertainties in the form factors are estimated to be ~8%. A further reduction in systematic uncertainties would 
require two-loop or nonperturbative current matching, finer lattice spacings, and light-quark masses at or very close 
to the physical values (which would also require increased lattice sizes). However, already at our current level of 
uncertainty in the form factors, the precision in phenomenological applications is primarily limited by the use of 
leading-order HQET for the heavy quark. 

To compare our form factor results to the literature, let us first consider the zero recoil point (_Ea = Wa, or 
(f — (Zmax)- There we obtain (matched to the MS scheme at ^ = mi,) 

Fi(gLx) = l-26(4)(10), 
J^2(gLx) - -0.288(25) (23), 
^^2(dax)M(9Lx) = -0.229(21)(11), (64) 

where the first uncertainty is statistical/fitting, and the second uncertainty is systematic. The ratio F2/F1 has 
previously been estimated by the CLEO collaboration using experimental data for the semileptonic Ac — > A e^Vg 
decay, assuming the same shape for Fi and F2 and ignoring Aqcd/^c corrections, to be F2/F1 = — 0.31(5)(4) 
This is consistent with our results, given the expected size of Aqcd/wc corrections, but significantly less precise. In 
Ref. [55], an earlier CLEO extraction of F2/F1 [53] was combined with the MIT bag model to obtain Fi(gmax) = 1-02, 
-^2('?niax) — —0.34. The authors of Ref. combined the CLEO data from Ref. [H] with the measured Ac lifetime 
to get i^i((j'max) — 1-21, -P2(9max) = —0.30, which happen to be quite close to our results. The sum rule calculation 
of Ref. [25j gave -Fi(9max) ~ 0.81, -F2(<Zmax) ~ —0.34, and ^2('?max)/^i('?max) ~ —0.42, in marked disagreement with 



(164|). 

At large recoil, leading-order soft-collinear effective theory predicts that Fi becomes equal to a soft form factor ^a, 
while F2 vanishes [27J [36l [37] . Using light-cone sum rules, the authors of Ref. [36] obtained ^a(9^ = 0) ~ 0.38. This 
is rather close to our results for Fi at that point [q^ = corresponds to E\ — ttia ~ 1.8 GeV), as can be seen in 



Fig. ( 14). However, we stress that our results are not reliable at such large values of i?A — tja, where we do not have 
lattice data and rely on extrapolation. 

Our form factor results can be used to make theoretical predictions for several observables in Af, -> A decays. 
As a first example, here we have calculated the differential branching fractions of these decays in the standard model. 
For i = fx, experimental data are already available from CDF 67J, and our calculation agrees with the data within 



uncertainties (see Fig. 17). The theoretical uncertainties in the branching fraction are dominated by errors of order 
AqcD/TOb and |p'|/mb associated with the use of leading-order HQET for the heavy quark, and by missing long- 
distance contributions to the decay amplitudes. At high recoil, there is also an unknown uncertainty associated with 
the extrapolation of the form factors. The long-distance contributions have not been calculated in lattice QCD, and 
need to be estimated using other approaches. The HQET uncertainties could be reduced by using higher-order lattice 
HQET, lattice nonrelativistic QCD, or a relativistic action for the b quarks. In such calculations one would have to 
deal with more complicated current matching and the full set of ten A^ — > A form factors. 

Another possible application of our form factor results to the phenomenology of Af, — A i^£~ decays is to study 
various angular distributions which depend on the baryon polarization and probe the helicity structure of the effective 
weak Hamiltonian. In particular, one should compute the angular distribution of the two-stage decay A;, — (A — > 
P'!T~) t^t~ for partially polarized A^, baryons. To make numerical predictions of these angular distributions for the 
LHC, the polarization of the A;, baryons produced in proton-proton collisions needs to be determined, for example 
using the method discussed in Ref. [3^. It remains to be seen whether the decay A^ — t' (A — >■ jjtt") will be 

competitive with B ^ K*{-^ KTr)i^£~ [64] in constraining new-physics models. 
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